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Abstract. In this paper we study cavity optomechanical systems in which the position of 
a mechanical oscillator modulates both the resonance frequency (dispersive coupling) and 
the linewidth (dissipative coupling) of a cavity mode. Using a quantum noise approach 
we calculate the optical damping and the optically-induced frequency shift. We find that 
dissipatively coupled systems feature two parameter regions providing amplification and 
two parameter regions providing cooling. To investigate the strong-coupling regime, we 
solve the linearized equations of motion exactly and calculate the mechanical and optical 
spectra. In addition to signatures of normal-mode splitting that are similar to the case of 
purely dispersive coupling, the spectra contain a striking feature that we trace back to the 
Fano line shape of the force spectrum. Finally, we show that purely dissipative coupling 
can lead to optomechanically-induced transparency which will provide an experimentally 
convenient way to observe normal-mode splitting. 
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1. Introduction 

Cavity optomechanical systems have a wide range of possible applications in precision 
measurement, quantum information, and fundamental tests of quantum mechanics [1-4]. 
In most optomechanical setups the coupling between the optical and mechanical degrees 
of freedom arises due to a displacement-dependent cavity frequency (dispersive coupling). 
Driving such systems at a frequency that is red-detuned from the cavity resonance can lead 
to cooling. This has been theoretically analyzed in Refs. [5,6], and the quantum ground- 
state has now been reached in several experiments [7, 8]. In the strong-coupling regime 
the optical and mechanical degrees of freedom hybridize. Normal-mode splitting has been 
predicted [9], and it was subsequently observed [10] in the optical output spectrum. If the 
system is probed with an additional probe field, the existence of the two normal modes can 
lead to destructive interference and a narrow transparency window at the cavity frequency 
[11]. For dispersive coupling this has been demonstrated experimentally [12, 13]. 

Recently, a different kind of optomechanics has been proposed [14]: a displacement- 
dependent cavity linewidth leads to a dissipative coupling between the mechanical and the 
optical degrees of freedom. Experimental realizations of this idea have been proposed in 
the microwave domain for superconducting resonators [14] and in the optical domain for a 
Michelson-Sagnac interferometer containing a moving membrane [15]. The ratio between 
dispersive and dissipative coupling is determined by the position of the membrane. A first 
experiment that demonstrated a dissipative coupling has been carried out with a microdisk 
resonator coupled to a nanomechanical waveguide [16]. It has been pointed out early on 
that dissipative coupling enables ground- state cooling outside the resolved- sideband limit 
and has potential applications in quantum-limited position measurements [14]. Moreover, 
squeezing of the mechanical state [17] and normal-mode splitting in response to a weak 
probe field [18] have been discussed in the context of the experimental setup of Ref. [16]. 
However, up to date, many properties of dissipatively coupled systems remain unknown. 

In this paper we study the general case of a cavity optomechanical system with both 
dispersive and dissipative coupling. After introducing our model in Section 2, we examine 
its mechanical and optical spectra in Section 3 and 4, respectively. For weak coupling we 
employ a quantum noise approach and calculate optically-induced damping and frequency 
shift of the mechanical oscillator, known as backaction damping and optical spring effect. 
In contrast to dispersive coupling, we find that dissipatively coupled systems feature two 
parameter regions of amplification and two parameter regions of cooling. We then present 
the exact solution to the linearized equations of motion for the general case of dispersive 
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as well as dissipative coupling. If the drive is red-detuned from the cavity resonance by 
the mechanical frequency, we find signatures of normal-mode splitting in the mechanical 
and the optical spectra. In the case of dissipative coupling a second feature appears which 
can be traced back to the Fano line shape in the force spectrum. In Section 5 we discuss 
optomechanically-induced transparency and find that this interference phenomenon can be 
observed for purely dissipative coupling. This could be useful to measure the normal-mode 
splitting. 

2. Model 

We consider an optomechanical system that consists of a mechanical oscillator with 
resonance frequency co m and a cavity mode with resonance frequency tu c . Dispersive 
coupling corresponds to a shift of the cavity resonance frequency due to the motion of 
the mechanical oscillator; dissipative coupling leads to a shift of the cavity damping rate 
due to the mechanical motion. The Hamiltonian of a dispersively and dissipatively coupled 
system is given by [14] (h = 1) 

= ui c a^a + co m tfb — 

The first term describes the cavity mode, where a) (a) are bosonic creation (annihilation) 
operators. The second term describes the mechanical oscillator, where & (b) are bosonic 
creation (annihilation) operators. The cavity has a linewidth k, and the mechanical 
oscillator is damped at a rate 7. The damping due to the optical and mechanical bath 
is described by H K and H^, respectively. The third term describes the optomechanical 
interaction with dimensionless coupling strengths A (dispersive) and B (dissipative), 
which are defined as derivatives with respect to the oscillator position x and given by 
Ak = — xq and Bk = ^r-%o, respectively. Here, x = (2mu m )~ 1 ^ 2 denotes the 
size of the zero-point fluctuations and m the mass of the mechanical oscillator; b q (K) are 
bosonic creation (annihilation) operators describing the optical bath coupled to the cavity 
and p denotes the density of states of the optical bath, treated as a constant for the relevant 
frequencies. B = corresponds to the well-investigated case of purely dispersive coupling 
and A = to the case of purely dissipative coupling. 



Ana* a + ^^/^^ Yl {f h i ~ b l^) + b )+H K + # 7 .(1) 
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To derive the Heisenberg equations of motion, we adapt the input-output formalism 
[19] to dissipative coupling. This leads to the following expression 

!— K kB /f f f \ 

= \/na in + -a + -— o + b Ma, (2) 




2 2 2 

where d in is the optical input mode [20]. The input-output relation is given by [15] 

flin - flout = -\f^a - (b + b 1 ) a. (3) 

2x V / 

Note that the last term only contributes for nonzero dissipative coupling and introduces 
an explicit dependence on the mechanical displacement as well as a nonlinearity into the 
input-output relation. 

Using a = (a + d)e~ luJdt , b = b + c, a m = (a in + ii n )e~ luJdt , and Eq. (2), we obtain 
the linearized equations of motion in a frame rotating at the drive frequency Ud 

a = - (iu m + |) c-yfyfj-HAK (a*d + aS^—y/n (a*| in - ~i\ ( n *d + , (4) 

'■} f . -r x\ <; k f ~ x \ i .— / B x \ £ / ~ k~ ^-EMx 
d = z A + Ak— d — 1 + B— d-y/H H £ in + iAkcl Bd - iQ— —.(5) 



x J 2 \ x J \ 2 i / \ 2 2 J xq 

In this expression, A = — oj c is the detuning between drive and cavity frequency, 
Vt = —iy/Ka- m is the strength of the coherent laser drive, and x = x (c^ + c) is the 
displacement of the mechanical oscillator relative to its steady-state position x. The 
thermal noise influencing the mechanical oscillator is described by the noise operators 
f} and if. The bath coupled to the mechanical oscillator is assumed to be Markovian and at 
a temperature T associated with an equilibrium phonon number nth — [^pi^m/ksT) — 
l] -1 , i.e. (if '(uj)fi(u')) = 2ir5(u + u!')n th and (fj(u)ff (u')) = 2it5(uj + u')(n th + 1) where 
ks denotes Boltzmann's constant. The operators £ in and describe the noise induced by 
the optical bath which is assumed to be vacuum noise, i.e. (£in(u;)£; n (u/)) = 2tt5(uj + u/). 

In Eq. (5) dissipative coupling B leads to a change in the damping rate k, whereas 
dispersive coupling A leads to a change in the detuning A. These shifts can be determined 
from the steady-state solutions of the classical equations of motion 

= b = - (iu m + b + iAn \d\ 2 - i^- (Qa* + Q*d) , (6) 



= d = % I A + An— ) a - 1 + -— )iQ - I 1 + B— ) -a, (7) 

X J \ 2 X V x oJ 2 
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Figure 1. (a) Electromechanical implementation of a dispersively coupled system (B = 
0). The resonance frequency of the LC resonator depends on the displacement-dependent 
capacitance Cq(x). The static output capacitor C\ determines the coupling strength 
between resonator and feed line, (b) Electromechanical implementation of a dissipatively 
coupled system [14]. Compared to (a) only the role of the capacitors is interchanged. This 
leads to a displacement-dependent coupling of the circuit to the feed line. As the total 
capacitance and thus the resonance frequency is displacement-dependent, the dispersive 
coupling is nonzero also in this case (i.e. A ^ and B ^ 0). 



where X — Xq (b + b*). These equations give rise to a static bistability, even if purely 
dissipative coupling, i.e. A = 0, is considered. This will be discussed elsewhere. 

Figure 1 shows electromechanical implementations of a dispersively coupled system 
and a system that includes both types of coupling [14]. In the case of purely dispersive 
coupling the interaction of the optical bath and drive with the mechanical oscillator is 
mediated by the cavity. Dissipative coupling also leads to a cavity-mediated influence on 
the mechanics, but in contrast to dispersive coupling it is proportional to the drive strength 
f2 instead of being proportional to the intra-cavity amplitude a (c.f. Eq. (4)). In addition, 
for dissipative coupling the mechanical displacement enters directly in the coupling to the 
optical bath, see Fig. 1 (b), and thus the mechanical oscillator is also directly influenced 
by the optical bath. This will lead to several new features in dissipatively coupled systems. 

3. Mechanical spectrum 

In the following we investigate the properties of the coupled system by discussing some 
of its fluctuation spectra. These spectra have the form Sk q (co) = J dt(y{t)q{Q))e tult = 
J du' (co)q(uj')) /2n and, in our case, require the solutions of Eqs. (4) and (5). 
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In this section we will focus on the properties of the mechanical spectrum S cc (cu) 
considering both weak and strong coupling. Applying a quantum noise approach we 
calculate the modifications of the mechanical spectrum due to weak dissipative coupling. 
We derive the optical damping and the optically-induced frequency shift from the weak- 
coupling force spectrum and compare them to the results in case of purely dispersive 
coupling. To go beyond weak coupling, we then present the exact solutions of the 
linearized equations of motion for the general case of dispersive and dissipative coupling. 
Finally, we calculate the mechanical spectrum and find signatures of normal-mode 
splitting and an additional feature originating from the modified force spectrum. 

3.1. Quantum noise approach 

Although the linearized equations of motion can be solved exactly, c.f. Eqs. (10) and (11), 
additional insight can be gained using a weak-coupling approach. In this subsection we 
discuss how the mechanical spectrum S cc (tu) is modified if only a small optomechanical 
coupling is considered. Recall that in absence of optomechanical coupling the mechanical 
oscillator is still coupled to the mechanical bath. Thus the oscillator is damped at a rate 
7 which leads to a mean phonon number in thermal equilibrium, n th . Together with 
the resonance frequency co rn , these quantities determine the mechanical spectrum S cc (co) 
which is a Lorentzian of width 7 (FWHM) with a peak at — tu m and an area of 27rn t h. 

For small coupling we treat the influence on the mechanical oscillator as a quantum 
noise source, inducing transition rates between neighbouring phonon number states, 
r n -m+i an d r n _>. n _i, given by Fermi's Golden Rule. Defining an amplification rate 
T t = T n ^ n+ i/{n + 1) = xlS FF (-uj m ) as well as a cooling rate T l = r„^ n _i/n = 
^Ff( u m) independent of the phonon number, both and are determined by the 
weak-coupling force spectrum Sff(u) [20]. These rates lead to an optically-induced 
damping 7 opt = — and a minimal phonon number n opt = Sff(— w m )/7o P t- In 
the presence of the mechanical bath and the optomechanical coupling this results in a total 
damping 7 to t = 7 + 7o P t which determines the new width of the Lorentzian describing the 
mechanical spectrum S cc (tu). Furthermore, the additional damping leads to a steady-state 
mean phonon number n osc = (771th + 7 op t^o P t) / (7o P t + 7)> thus the area of the Lorentzian 
is changed. Finally, optical damping affects the effective spring constant corresponding 
to a shift of the mechanical frequency given by 5cu m = J duSpF^)^-! '(w m — cu) — 
l/(co m + uj)]/2it [5]. Thus, the modifications of the mechanical spectrum due to weak 
optomechanical coupling can be described by the parameters % pt , n osc and 5u m . 
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So far our considerations do not explicitly depend on the type of the coupling. 
However, the force spectrum S FF (u) contains this information, i.e. its shape depends on 
the applied coupling. To calculate S FF (u) we use the backaction force operator [14] 



B * * B 

Fxq = AK{a*d + ad)) + i— — aCi n ) — + 



(8) 



which can be determined from the interaction part of the Hamiltonian, i.e. the third term 
of Eq. (1). Substituting Cl = —ia(iA — k/2) derived from the steady-state solution of the 
uncoupled Eq. (7), the force spectrum is given by [14] 

+oo 



S ff (uj) 



dJ 

27 




1 + -iA 




Xc{UJ) 



(9) 



IXc(w)| ! 



2A 



2Ak 
B 



In the general case of dispersive and dissipative coupling the result is a Fano line shape 
which reduces to a Lorentzian in absence of dissipative coupling, i.e. B = 0. As discussed 
in Ref. [14] the Fano line shape originates from an interference effect between the two 
ways of interaction with the mechanics. These two processes act as two noise sources 
influencing the mechanical oscillator, and lead to the two terms inside the absolute value 
in Eq. (9): The constant first term accounts for the direct interaction between optical bath 
and mechanical oscillator and represents coupling to a continuum. In contrast, the second 
term is filtered by the cavity response x c (w) and arises due to the influence of the cavity. 
The interference of these two contributions, the direct action of the optical bath and its 
cavity-mediated influence, gives rise to the Fano line shape. Purely dispersive coupling 
leads only to a filtered, cavity-mediated influence, i.e. the mechanical oscillator is only 
influenced by a single optical noise source, and no interference can occur. 

Figure 2 (a) shows the optical damping for purely dispersive and purely dissipative 
coupling. Since B = means that the force spectrum is a Lorentzian, the optical damping 
7o P t is given by the difference of two Lorentzians. Choosing A w — cu rn maximizes 
the optical damping rate. In contrast, since the force spectrum S FF (tu) is a Fano line 
shape for dissipative coupling (purely or in addition to dispersive coupling), the optical 
damping rate is modified [14]. The maximum is shifted farther away from the mechanical 
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Figure 2. Optical damping 7 opt (a) and optically-induced frequency shift 5ui m (b) as a 
function of detuning A. The dashed, green lines show the result for purely dispersive 
coupling (Aa = 0.4, B — 0), the solid, black lines show purely dissipative coupling 
(A = 0, Ba = 0.4). Blue (red) areas in (a) indicate cooling (amplification). The sideband 
parameter is w to /k = 3. 



resonance frequency u m and for |A| ^> k the optical damping rate decreases more slowly 
than a Lorentzian and is proportional — 1/A. Furthermore, for typical parameters we find 
two regions where the optical damping is positive, thus providing cooling, as well as two 
regions with negative 7 opt , leading to instability if the total damping 7 tot = 7 + 7 opt < 0. 

Figure 2 (b) shows the optically-induced frequency shift 5co m for purely dispersive 
and purely dissipative coupling. Dispersive coupling and cooling at A = — ui m allows 
for a vanishing frequency shift 5u m = 0. In contrast, dissipative coupling leads to a 
nonzero frequency shift §tu m at A = — u m . It remains small for small detunings only, 
for large values of A it increases linearly. Note that this linear dependence is due to the 
fact that we fix the number of photons inside the cavity \d\ 2 , which implies that the drive 
strength f2 has to increase with the detuning A. Since dissipative coupling has a component 
proportional to f2, the effective dissipative coupling strength is increased. Fixing the laser 
power instead, the intra-cavity amplitude a decreases as 1/A and 5u> m ~ (5|a|) 2 A/2 goes 
to zero in the limit of large detunings |A| ^> uj m . 

One possible choice to achieve cooling with dissipative coupling is A « — co m . Figure 
3 shows that, for both purely dispersive and purely dissipative coupling, A w — to m leads 
to a strong decrease of the phonon number n osc . This is not surprising since the optical 
damping 7 opt is maximized close to this detuning in both cases. However, large optical 
damping alone is not sufficient to achieve the best cooling results in the sense of smallest 
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Figure 3. Mean phonon number n osc as a function of detuning A in case of (a) purely 
dispersive coupling and (b) purely dissipative coupling. The solid lines show the result for 

(a) Ad = 0.01 and (b) Ba = 0.01, respectively, the dashed lines show (a) Ad — 0.05 and 

(b) Ba = 0.05, and the dot-dashed lines show (a) Aa = 0.3 and (b) Bd = 0.3. Other 
parameters are w m /« = 3, cj m /j — 10 5 , and n t jj = 100. Hatched areas indicate unstable 
regions due to the criterion 7 tot < 0. 



n osc . Notably, dispersive coupling at this detuning leads to smaller n osc despite the larger 
optical damping rate of dissipative coupling, except for very small coupling strengths. 
This is due to a larger n opt which also contributes to the mean phonon number n osc . Thus 
it is of particular interest to achieve n opt as small as possible. For dispersive coupling this 
is linked to reaching the resolved- sideband limit [5]. In presence of dissipative coupling 
the Fano line shape of Sff(u) leads to an optimal detuning A opt = u rn /2 + kA/ B where 
^o P t = [14]. More generally, the Fano line shape of Sff{u) allows to find a detuning 
A (w) = — u/2 + kA/B for each uj, such that Sff{u) = 0, c.f. Eq. (9). Hence, A opt = 
A (— u! m ) implies S F f(— w m ) — and this in turn implies n opt = S F f(— w m )/7 op t = 
and = 0. Therefore the optomechanical coupling induces a cooling rate but no 
amplification rate and ground-state cooling can be achieved if the drive strength is 
sufficiently large or the intrinsic damping 7 small enough. Fortunately these conditions are 
independent of the sideband parameter co m / k and ground- state cooling can be performed 
in the unresolved- sideband regime that is easier to reach experimentally. Finally, since the 
optimal detuning is part of the second cooling region, A = A opt is far from maximizing the 
optical damping rate, see Fig. 2 (a). Compared to the values of 7 opt achieved at A = — cu m 
for either dispersive or dissipative coupling, the optical damping rate at A = A opt is rather 
small. Therefore, to achieve considerable cooling despite the poor cooling rate, stronger 
coupling or smaller intrinsic mechanical damping 7 is required. Figure 3 (b) shows that 
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the advantage of n opt = combined with sufficiently strong drive, allows to exceed the 
cooling result of dissipative coupling at A = — cu m . 

Note that in the case of purely dissipative coupling, i.e. A — 0, the optimal detuning 
A opt = to m /2 corresponds to a blue detuned drive laser. On the contrary, driving a 
dispersively coupled system (B = 0) with this detuning would lead to amplification rather 
than cooling. 



3.2. Strong coupling 



In this subsection we want to investigate the mechanical spectrum beyond weak coupling. 
Thus, we first calculate the exact solutions of Eqs. (4) and (5) for the general case of 
dispersive and dissipative coupling. Then, using these results, we derive the mechanical 
spectrum S cc (cu) and discuss its features, i.e. compare the case of purely dissipative 
coupling to the case of purely dispersive coupling. 

To proceed from Eqs. (4) and (5), we solve the uncoupled (A = B = 0) classical 
equations (6) and (7) in steady-state and find x = and = — id (iA — |) . Using these 
results and solving the coupled equations of motion (4) and (5) in Fourier space, we obtain 



} 



(10) 



Xm 1 (~ UJ ) a*a(u)£ in (uj) - aa*(-uj)H(uj) 



B 



[4 1 (-wWH + ^*)] (id 



2iu m 



B 



— a[uo\ 



\a\ 2 a(uj)£ in (uj) - a 2 a*(-uj)i} n (uj) 



Here we have used the cavity response function Xc(u) = [ K /2 — i{w + A)] 1 , the 
response of the mechanical oscillator Xm{u) — [7/2 — i{u — oo m and Jf{u) = 
-uj) + 2cu m E(a;). Furthermore, we have defined the optomechanical self- 



*-!/ 



energy S(w) and a(co) as 



(12) 
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where 



s aM = -i(An\a\) 2 [xc(w) - Xt(-u)} 

= i (f ) 2 |a| 2 [ Xc ( w ) (zA + f) 2 - xt(-co) (<A - f) ; 
E^w) = 5i K |a| 2 [ Xc (a;) (zA + f ) - ^(^w) (zA - f )] 



(13) 



and 



(14) 



ag(u) = f - f xM (iA+f). 

Note that the Fourier transformation was applied such that <^(w) = [Q{~ for all 
operators. In the purely dispersive case (_B = 0) the above definition of the optomechanical 
self-energy E(u;) reproduces the notation used in [5]. Defining Eg(a>) and E^^(o;) in a 
similar fashion, we can deduce the optical damping 7 opt = — 2Im[E(u; m )] and frequency 
shift 5tu m = Re[E(u; TO )]. Note that this means that E^(o;) differs from the definition 
in [14] by a factor of — 2iuj m x* m {— w). 

Remarkably, the exact solutions of the linearized equations of motion, Eqs. (10) 
and (11), have the same structure for both types of coupling. Apart from an additional 
contribution proportional to B in Eq. (11), differences are hidden in the functions E(o>) 
and a(u). The additional dependence of E_g(u;), E^g(u;) and a^(u) on the detuning A 
arises since, for dissipative coupling, the equations of motion (4) and (5) contain a term 
proportional to the drive strength f2. The constant term in a§(u) is due to the direct 
interaction between the optical bath and the mechanical mode. 

Finally, the fluctuations of the optical output are obtained by using a out = (a out + 
iout)e~ tuJdt and linearizing the input-output relation (3), £ in — £ out = —^/nd— y/K,Bdx/2x . 
Then, with Eqs. (10) and (11), we find 

+ [1 - kxM] 6» + 2%K ^^ ] [\a\ 2 a(u)U") - a 2 a*(-u)il(u) 
Using Eq. (10) we can calculate the mechanical spectrum 



+oo 



— (tf{u})c{u)) = r / 

^7T \jV{u)\ 



Scc{u)= I —{c j {uj)c{uj)} = , — 2 , (16) 



where a th (uj) = |£(w)| 2 (n th + 1) + \Xm( u ) + i£(w)| 2 n th and a opt (u) = 

IXm 1 ^)] 2 S ff {uj)x1/ k. This result is valid for purely dispersive, 
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Figure 4. (a) Logarithm of the mechanical spectrum S cc (u)k as a function of detuning 
A. Parameters are u m /n = 3, ui m H = 10 5 ' n th = 100, A = 0, and Ba — 0.4. The 
dark regions indicate regions of instability obtained from a numerical calculation. The 
green curve gives half of the total damping rate, 7 to t /2, obtained from the quantum noise 
approach with the origin shifted to (—2.5,0). The dashed lines show the real part of the 
eigenvalues of the Hamiltonian (17). (b) S cc (uj)k for detunings A/uj m = 0.55, 0.5, 0.45 
(from top to bottom), (c) S cc (u>)k for detunings A/u) m = —0.9, —1, —1.1, —1-2, —1.3 
(from top to bottom). 

purely dissipative and both types of coupling but has the same form as found in the case 
of dispersive coupling only [5]. For B = the result coincides with [5]; setting A = the 
result coincides with [14]. 

Figure 4 (a) shows the mechanical spectrum for strong dissipative coupling. Dark 
areas indicate regions where the solutions of the linearized equations of motion are 
unstable. This was numerically tested for the parameters used in Fig. 4 and coincides 
with the regions where the total damping rate 7 to t from the weak-coupling approach is 
negative. Whereas dispersive coupling leads to one unstable region for blue detuning, 
dissipative coupling can lead to a second unstable region for red detuning in addition to 
an unstable region for blue detuning. A third unstable region exists for even stronger drive 
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or large red detuning. This is not predicted by the behaviour of the optical damping rate, 
i.e. it appears although 7 tot > 0. 

Focusing on the stable regions, we find two prominent features. First, at A = A opt = 
co m /2 a strong decrease of the phonon number (n) = / du:S cc (uj) /2ir can be observed. As 
mentioned in Subsection 3.1, this detuning is associated with cooling [14] and a special 
case of the strong modifications of the mechanical spectrum at A (u>). If A = A (w), the 
force spectrum Sff(u) vanishes, which means that at this frequency cu only the first term 
of Eq. (16), cr t h(w), contributes to S cc (ui). Furthermore, Fig. 4 (b) shows that, apart from 
the main peak close to the mechanical resonance to = —co m , there is a broad contribution 
at a second frequency arising from a opt (<jj) /\<yY (u)\ 2 . It is given as a trade-off between 
the maximum of the Fano line shape of the mechanical spectrum Sff{u) at to = ~ 4 ^ +k2 
and the peak of \ jV{uj)\~ 2 at to = —cu m . It is this contribution, away from the mechanical 
resonance frequency, that finally limits the cooling due to its increasing relevance with 
increasing coupling strength. 

The second feature is found at A = — uj m . Similar to the case of dispersive coupling, 
we find normal-mode splitting even though slight quantitative differences appear. In the 
following we use a simplified Hamiltonian to find an approximation that describes the 
splitting. Recall that dissipative coupling leads to two terms in the equations of motion (4) 
and (5). We neglect the term proportional to the damping rate, i.e. the direct influence of 
the optical bath on the mechanical oscillator, and only take the effect proportional to the 
drive into account. Furthermore, we use the rotating wave approximation and neglect the 
fast rotating terms d)$ and dc. Then, in the rotating frame, the simplified, non-hermitian 
Hamiltonian is given by 



H = - (A + i^j d ] d + (u m - C^j c ] c + 



B — - Aan ) cd) + H. 



(17) 



Note that using this approximation, the difference between purely dispersive and purely 
dissipative coupling only depends on whether the drive strength or the intra-cavity 
amplitude a is fixed. Fixing Cl for purely dissipative and a for purely dispersive coupling 
leads to similar results. Instead fixing one parameter for both types of coupling, as done 
here with a variable drive strength Vt and a fixed a, leads to modifications of the splitting 
due to an additional dependence on the detuning A. Since it is not possible to fix both 
Vt and a at the same time, mixed coupling will always lead to A-dependent modifications 
arising from either the dissipative or the dispersive term. 

In the general case of dispersive and dissipative coupling, the eigenvalues of the 
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simplified Hamiltonian (17) can be calculated as 

p _ . 7 + k oj m — A 



IQ(A 2 k 2 - ABAk) + B 2 (4A 2 + k 2 ) . 

The energies corresponding to the two modes are the real parts of these eigenvalues E±, 
whereas the imaginary parts contain information about the associated linewidths. We show 
the real parts of the eigenvalues (calculated for A = 0) in Fig. 4 (a) and, despite the 
simplifications, the energies fit the peak position of the spectrum very well. Differences to 
purely dispersive coupling arise since the dispersive coupling matrix element is constant 
for fixed values of the cavity amplitude a. In contrast, the dissipative coupling matrix 
element depends on the drive strength fi, which is a function of detuning A if a is 
fixed. This affects the curvature of the modes and leads to A-dependent width of 
the splitting. Moreover, the eigenvalues of the simplified Hamiltonian indicate that, in 
case of purely dissipative coupling, the splitting is no longer minimal at A = — tu m . 
Neglecting the damping terms in the Hamiltonian (17) the minimal splitting occurs at 
A = — tu m /(l + B 2 \a\ 2 ). Figure4(c) shows in detail how the single peak at the mechanical 
frequency is split due to the optomechanical coupling. 

We further investigate the eigenvalues E± from Eq. (18) to clarify at which coupling 
strength normal-mode splitting appears in the mechanical spectrum S cc (tu). Figure 5 (a) 
shows that Re[E7±] coincide well with the peak positions of the mechanical spectrum as a 
function of coupling strength Bd. If A = — u m , small coupling corresponds to degenerate 
energies, i.e. Re[£L] = Re [.EL]. In this case, the argument of the square root in Eq. (18) is 
real and negative. Thus, for small coupling, the root contributes only to the imaginary part 
of the eigenvalues and affects the linewidths given by k and 7 respectively. With increasing 
coupling strength the linewidths approach their mean value (k + j)/2, which is reached 
where the root becomes zero. Then the modes Re[.EL_] start to split whereas the linewidths 
remain unchanged. In case of purely dissipative coupling at A = — co m , the critical 
coupling strength where mode-splitting starts is given by B\a\ = (k — 7)/^/4a;^ + k 2 . 

Figure 5 (a) shows normal-mode splitting in case of purely dispersive coupling. 
For the set of parameters used, it starts at a larger coupling strength than the splitting 
obtained for purely dissipative coupling. Note, however, that this depends on the sideband 
parameter co m /K, since the critical dispersive coupling strength at A = — tu m is given 
by A\a\ = (k — 7)/(4ac). Thus, for w^/zt 2 < 15/4, dispersive coupling would lead to 
normal-mode splitting at a smaller coupling strength than dissipative coupling. 
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Figure 5. (a) Real part (black curves) and imaginary part (green curves) of the eigenvalues 
E± calculated from the Hamiltonian (17) as a function of the coupling strength. Solid 
(dashed) lines indicate purely dissipative (dispersive) coupling. The plot is underlaid with 
the mechanical spectrum S cc (uj)k as a function of coupling strength Ba for A = — oj m and 
A = 0. Other parameters are uj m j 'k = 3, u! m /j — 10 5 , and n t h = 100. (b) Mechanical 
spectrum S cc (u)k at A = — oj m for different coupling strengths Ba as well as A = 0. 



Finally note that if A ^ —u> m , the root in Eq. (18) is complex valued and the modes 
start with a finite energy separation from the uncoupled case. 

4. Optical output spectrum 

The optical spectra, especially the optical output spectrum, are experimentally easier 
accessible than the mechanical spectrum. Thus we use the full solutions (11) and (15) 
to calculate the cavity and the optical output spectrum. Purely dispersive coupling allows 
interaction between the mechanical element and the optical output only via the cavity, such 
that (w) = nSdd(u). Note that this is no longer the case for dissipative coupling since 
there is direct influence of the mechanical oscillator on the output which is not mediated 
by the cavity. We find 



Strong-coupling effects in dissipatively coupled optomechanical systems 



16 



and 





(20) 



= K 



a(-u)-B/2\ 2 



Apart from the factor k, these two spectra differ by the subtraction of the constant term 
from a(—u}). Recalling the definition of a(u), we find that this means, that dissipative 
coupling contributes to the cavity spectrum only at frequencies filtered by the cavity 
response \Xc{~ w )| 2 - This leads to the enhancement of the lower sideband for A < 
and of the upper sideband if A > 0, similar to the case of dispersive coupling. Due to 
the direct influence of the mechanical oscillator on the optical output, dissipative coupling 
leads to a contribution to the output spectrum S^f (w) that is not filtered by the cavity 
response. This is hidden in the definition of a(—u) in Eq. (14). 

The optical output spectrum is connected to the displacement spectrum S xx (oS) via 
Sddi 03 ) = S ff (—uj)S xx (uj) [14]. Thus, it is possible to observe the features of the 
mechanical spectrum S cc (lo)k in the optical output spectrum. As shown in Fig. 6 (b) 
we recover normal-mode splitting at A = —co rn and find modifications of the optical 
output spectrum S^(oj) for A = A (±w). First, we can see the influence of A on the 
mechanical spectrum S cc (co) at A = A (+w). Moreover, there is also the direct influence 
through the weak-coupling force spectrum Sff(—u), i.e. the optical output spectrum 
becomes exactly zero if A = A (— u). Figure 6 (a) shows S^f (w) in case of purely 
dispersive coupling for comparison. Normal-mode splitting can be observed as well, but 
for purely dispersive coupling the detuning A has no special role. Note also the different 
instability regions of the optical output spectrum, depending on the type of coupling. 

For mixed coupling (i.e. A ^ and B ^ 0) the new features of dissipative coupling 
are modified but do not disappear. In particular, there is a detuning A (u/) such that 
Sff(co) = 0. However, its offset kA/B depends on the ratio of the couplings and leads 
to a shift of A compared to the purely dissipative case. Furthermore, mixed coupling 
modifies the regions where 7 tot < in the weak-coupling approach and the corresponding 
changes of the unstable regions are captured by the numerical calculation as well. 

5. Optomechanically-induced transparency (OMIT) 

In this section we will investigate the response of the optomechanically-coupled system 
to a weak probe field, and we show that purely dissipative coupling, i.e. A = 0, leads to 
optomechanically-induced transparency. This is also a convenient way to observe normal- 
mode splitting (NMS). We compare our findings to the purely dispersive case and give 
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a>/u) m co/co, 



Figure 6. Logarithm of the optical output spectrum Sffi (w) for (a) purely dispersive 
coupling (Aa — 0.4, B = 0) and (b) purely dissipative coupling (A = 0, Ba — 0.4). 
Other parameters are the same as in Fig. 4. Dark regions indicate regions of instability. 
The white line indicates an exact zero of S^fa). The dashed line indicates where the 
mechanical oscillator experiences dissipative cooling associated with the white line in 
Fig. 4 (a). 

an appropriate approximation that holds in the general case of dispersive and dissipative 
coupling, i.e. A ^ and B ^ 0, in the resolved- sideband regime. 

The probe field of frequency co p is assumed to be weak compared to the drive field, 
i.e. its optomechanical coupling can be neglected. Thus it is sufficient to account for 
the probe laser by changing the optical input mode £in(w) in an appropriate way and 
neglecting additional coupling terms. In Section 3 and 4 the operator £ in denoted vacuum 
fluctuations only, now it contains the probe field such that £i n (t) = £, vac (t) + d p e~ lSt with 
(lin(*)) — d P e~ l5t . Here <f vac describes the vacuum fluctuations of the optical bath and 
5 = Up — Lu d denotes the detuning between probe and drive laser. 

We investigate the response to the probe field by evaluating the expectation value of 
the optical output mode (15). Since (17) = (tf) = 0, the result is of the form 

/+°° A, . 
— {U t (u))e-^ = A-e~ m + A+e t5t . (21) 

Recall that all calculations are done in a frame rotating with — w d , thus the optical output 
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contains terms rotating at three frequencies: — Ud (drive frequency), —5 — uJd = —lo p (anti- 
Stokes field) and 5 — Ud = co p — 2ud (Stokes field). The contribution at the drive frequency 
is not contained in Eq. (21) since we treated the coherent part of the drive separately with 
Eqs. (6) and (7), i.e. £ out only describes the fluctuations around the strong drive field. A" 
and A + are the complex amplitudes of the anti-Stokes and Stokes field and are given by 

\a\ 2 a(Sy 



A' 



1 — KXc(8) + 2lKUJ ri 



(22) 



A + = -2iKW m iri, X robe- (23) 



a 2 a*(5)a(- 5) J# 

The anti-Stokes field rotates with the probe frequency — tu p , thus A~ is the amplitude of the 
original probe field modified due to interference with anti-Stokes scattered light (5 > 0) 
from the drive field. Furthermore, A + is the amplitude of the output field component 
rotating at a frequency cu p — 2tUd that is created by the optomechanical coupling, i.e. Stokes 
scattering (5 > 0) of drive photons. 

Focusing on the anti-Stokes contribution at A = — cu m where NMS appears in the 
mechanical spectrum S cc (co), Eq. (22) consists of three contributions to the amplitude A~: 
The constant first term accounts for the initial probe field. The second term represents 
the influence of the uncoupled cavity. Finally, the third term is nonzero only for nonzero 
coupling and contains the influence of both dispersive and dissipative coupling. 

Using homodyne detection, different quadratures of the anti-Stokes field can be 
investigated experimentally. Figure 7 shows the real part of the anti-Stokes amplitude, 
Re[v4~]. In absence of optomechanical coupling the cavity leads to a Lorentzian- shaped 
dip of width k associated with the second term in Eq. (22). For nonzero coupling the 
third term in Eq. (22) modifies the amplitude due to scattering processes from the drive 
to this frequency. These processes are suppressed away from the mechanical resonance, 
thus striking modifications occur only for S ~ ±w m . There, an upper or lower sideband 
— Ud ± co m is created and its frequency coincides with the probe frequency — cu p , which 
gives rise to interference effects [12, 13]. 

At S = —u rn , scattering from the drive laser is not suppressed by the mechanical 
response, but in the resolved sideband regime, i.e. co m k, this process is highly off- 
resonant with respect to the cavity frequency. Thus, the effect at this frequency is small. 
As shown in the insets of Fig. 7, dissipative coupling leads to a larger contribution at 
this detuning 5 ~ — cu rn than purely dispersive coupling. This originates from the direct 
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Figure 7. The real part Re[A~ / d plo b e ] of the response at the probe frequency — uj p as a 
function of the detuning between probe and drive field S for A = — uj m . (a) and (b) show 
the case of purely dispersive coupling with Ad = 0.1 and Ad = 0.4. (c) and (d) show 
purely dissipative coupling with Bd = 0.1 and Bd = 0.4. Other parameters are the same 
as in Fig. 4. The green dashed line shows the result in absence of coupling (A = B = 0). 
The insets show a magnification around S = —u m . 

interaction between optical bath and mechanical oscillator: It gives rise to a constant 
contribution to a(S) 2 in Eq. (22), i.e. a term not filtered by the cavity response function. 

In contrast, if 5 = +co m and A = —co m , the probe frequency co p coincides with the 
cavity resonance co c , giving rise to more prominent effects. The optomechanical coupling 
leads to a narrow peak enclosed by a broad dip that appears also in absence of coupling. 
For small coupling as shown in Figs. 7 (a) and (c), the width of this peak is given by the 
width of the mechanical resonance. The mechanical linewidth, in turn, is given by the 
intrinsic damping 7 and broadened with increasing coupling strength due to the additional 
optical damping 7 opt . In the case of sufficiently strong coupling, see Figs. 7 (b) and (d), 
the two modes are separated by a peak that has a width comparable to or larger than the 
width of each of the modes. The splitting increases for stronger coupling. 

This general behaviour is shared by dissipatively and dispersively coupled systems, 
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but there are small differences: First, the width of the splitting in case of purely dissipative 
and purely dispersive coupling depends differently on the respective coupling strength. 
Second, there is an increasing asymmetry between the two modes in case of dissipative 
coupling, whereas purely dispersive coupling leads to a splitting into two anti-peaks that 
remain similar over a larger range of coupling strengths. 

In analogy to the treatment of purely dispersive coupling [12], we assume that only 
anti-Stokes scattering occurs. This can be described by simplified equations of motion 
where we neglect coupling to d) and & in Eqs. (4) and (5). As a result (£ ut(£)) is still °f 
the form of Eq. (21), but with new coefficients A+ pprox = and 



a- 

approx 



\a\ 2 a{5) 2 



1 - KXc(<5) - k 



Wprobe 

(24) 



where £(u>) = + + ^ab( u )- H ere > = — i{AK\a\) 2 Xc{u), £g(u;) = 

i{B/2) 2 \a\ 2 Xc{u){i A + k/2) 2 , and t M (uo) = BAk \a\ 2 x c (w)(z'A+k/2) denote only the 
parts of the originally defined self-energies with weight at 5 ~ +u m . In the case of B = 
the approximation is the same as in Ref. [12]. For both dispersive and dissipative coupling, 
the approximation is valid in the resolved sideband regime, e.g. it does not reproduce the 
feature at S = —co m which becomes more important if uj m is of the order of k or the 
coupling becomes too strong. 



6. Conclusion 



We have presented a detailed study of optomechanical systems featuring both dissipative 
and dispersive coupling. For weak coupling we have employed a quantum noise approach 
to calculate the optical damping and the optically-induced frequency shift. Surprisingly, 
there are two regions leading to cooling and two regions leading to amplification. This 
is a consequence of the Fano line shape in the force spectrum which is absent for purely 
dispersive coupling. In the strong-coupling regime we have calculated the mechanical and 
the optical spectra from the exact solution to the linearized equations of motion. Similar to 
purely dispersive coupling, normal-mode splitting appears for sufficiently strong coupling. 
Nonzero dissipative coupling additionally leads to a striking feature which originates from 
quantum noise interference. Finally, we have found purely dissipative coupling can lead to 
optomechanically-induced transparency which will be an experimentally convenient way 
to observe normal-mode splitting. 
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